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ABSTRACT 

In this paper, we have pressed a method of evaluation of energy 
per particle (E/N) of two dimensional Bose Gas using three 
potential namely HD (Hard-disc potential), SD (Soft disc 
potential), PP (Zero-range pseudopotential). By using a zero-range 

pseudopotential we have estimated na^p = 0.o4 as the critical 

density at which the system becomes unstable against cluster 
formation. 
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1. INTRODUCTION 

Degenerate low-dimensional gases 
are presently attracting considerable interest 
as model system to investigate beyond 
mean-field effects and phenomena where 
many-body correlations and thermal and 
quantum fluctuation play a relevant role. 1 In 
two dimensions (2D), it is well known that 
thermal excitation destroy ling-range order 
and Bose-Einstein Condensate (BEC) cannot 
exist in large bosonic systems at finite 



temperature. Nevertheless, a defect mediated 
phase transition from a high temperature 
normal fluid to a low temperature superfluid 
was predicted by Berezinski, Kosterlitz and 
Thouless (BKT) 3 and was observed in thin 
films of liquid 4 He. 

At zero temperature, fluctuations are 
instead negligible and BEC is possible. The 
ground-state energy per particle of a 
homogeneous dilute Bose gas in 2D was 
first calculated by Schick 5 and is given by 
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Where m is the mass of the particle, n is the 
number density, and a 2 D is the 2D scattering 
length. The above result holds for small 

values of the gas parameter, na\ D « 1, and 

can be obtained within a mean-field 
approximation using the coupling constant 
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The presence at low energy of a 
density-dependent coupling constant is a 
peculiar feature of 2D systems. 

Bose gases in quasi-2D have been 
realized with spin-polarized hydrogen on a 
liquid-helium surface 6 and with alkali- metal 
atoms confined in highly anisotropic 
harmonoic traps. 79 In these systems, the 
transverse motion of the atoms is frozen to 
zero-point oscillations resulting in fully 2D 
kinematics, but the interatomic scattering 
processes are still 3D. In fact, the 3D s-wave 
scattering length a 3D is always much smaller 
that the characteristic length of the tight 
transverse confinement. Theoretical studies 
of quasi-2D trapped systems predict the 
existence of a low-temperature quasi 
condensate phase 10 ' 11 whose properties have 
been investigated using mean-field 
approaches. 12 The nature of the transition in 
confined systems and whether it belongs to 
the BEC or BKT universality class is still an 
open problem both experimentally and 
theoretically. 

An important theoretical prediction 
for quasi-2D Bose gases in harmonic traps is 



the renormalization of the coupling constant 
due to the tight confinement in one direction. 
To logarithmic accuracy the result is given 

by 11 ' 13 

2V2^ 2 1 



m 



- + 



i 3D 



= ln\\ln(0)al\ 



(3) 



Where a z = I mW z is the oscillator 

length, fixed by the frequency 0) z of the 

tight transverse confinement, and n(o) is the 
2D density in the center of the trap. If 
a 3D «a z one can neglect the logarithm term 
in Eq. (3) and the resulting coupling constant 
is determined by the 3D scattering length. 
On the contrary, If a 3D »a z ,g 2 D becomes 
independent of the value of a 3D and a 
regime of pure 2D scattering is achieved 
with a 2 D=a z . This regime is highly 
interesting and can be realized in trapped 
gases by using a Feshbach resonance 14 to 
achieve large value for a 3D . 



2. MATHEMATICAL FORMULAE 
USED IN THE EVALUATION 

We consider a homogeneous system 
of N spinless bosons 2D described by the 
many-body Hamiltonian 



(4) 



Where r = x i + y j denotes the 2D 
position vector of the ith article and V(r) is 
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the two-body interatomic potential. We are 
different choices for the interatomic 
potential V(r). 



(i) Hard-disk (HD) potential defined by 



(5) 



0(r>a 2D ) 

in which case the hard disc diameter 
corresponds to the scattering length a 2D . 



(ii) Soft-disk (SD) potential defined by 
\v 0 {r<R) 



V SD {r): 



0{r>R) 



(6) 



in which case the scattering length is given 
in terms of the modified Bessel function 
I 0 (x) and its derivative I' Q (x) 



l 2D 



■ Rexp 



Io(RK Q ) 



RK 0 I' 0 (RK 0 ) 



(7) 



where K% = mV 0 / h 2 with V 0 >0. For a 

finiteV ( , one has always R>a 2D; if V 0 — > +°°, 
the SD and HD potential coincide with 
a 2 D=R In the present calculations, the range 
of the SD potential is kept fixed to the value 
R=5a 2D and the height V 0 is determined 
through Eq. (7) to give the desired value of 

g2D. 

(iii) Zero-range pseudopotential (PP) 



defined by 

2xtf 8{r) 



V pp (r) = 



m In(qa 2D ) 



\-In(qr)r — 
or 



where the second term in the square brackets 
is a regularizing operator with q an arbitrary 
wave vector. The potential V pp (r) supports a 
two-body bound state with energy 

e b a- Ah 2 1 [ma\ D eY] 

And wave function 

f b {r) = K 0 (2r/{a 2D ey)), 

where K 0 (x) is the modified Bessel 
function and f = 0.577 is the Euler 
constant. 

The calculations with the repulsive 
potential V HD (r) and V SD (r) are carried out 
using the diffusion Monte Carlo (DMC) 
method. This technique solves the many- 
body time independent Schrodinger equation 
by evolving the function 

f{R,T) = ys T {R)y{R,T) 

in imaginary time T — it I h according to the 
time dependent Schrodinger equation 

-^^- = -DV 2 R f(R,T) + DV R [F(R)f(R,t)] + 



[E L (R)-E ref ]f(R,T) 



(9) 



Here ^(T?,?") denotes the wave function 
of the system and l/f T (R)is a trial function 

used for importance sampling used for 
importance sampling. In the above equation 

R = {r v ..,r N ),E L {R) = ¥T {R)- 1 H ¥T {R) 

Denotes the local energy, 

F{R) = 2 ¥t {R)- 1 V r¥t {R) 

is the quantum drift force, D -h 2 I (2m) 
plays the role of an effective diffusion 
constant, E ref is a reference energy 



(8) introduced to stabilize the numeric. The 
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energy and other observables of the ground 
state of the system are calculated from 
averages over the s\asymptotic distribution 
function >°°). Apart from 

statistical errors, the DMC method 
determines the ground-state energy of a 
system of N bosons exactly. 

In the case of the pseudo potential 
V pp (r) Eq.(8), the gas like state is not the 
ground state of the system, but is a highly 
excited state of the Hamiltonian, We except, 
however, that for small values of the gas 

parameter, na\ D « 1, the gas like state is 

stable against formation of many-body 
bound states. The calculations with the PP 
potential are performed at the level of the 
variational Monte Carlo (VMC) method to 
avoid the technical difficulties that arise in 
the study of excited states in the DMC 
method. Nevertheless, the accuracy achieved 
by using this variational approach is very 
high. The VMC technique is based on the 
variational principle 



\ H \Wt) 



(Yt\Yt) 



>E 



which is here applied to excited states of the 
Hamiltonian H. For a trial wave function 
y/ T with a given symmetry, the variational 

estimate provides an upper bound to the 
energy E of the lowest excited state of the 
Hamiltonian H with that symmetry. 

In all the calculations the trial wave 

function y/ T [R) is taken of the Jastrow 
form 



en) 



Where the correlation factor f(r) is chosen as 
the solution of the two-body Schrodinger 
equation subject to the boundry conditions 
f(r>d)=l and f'(r>d)=0 on the wave 
function and its first derivative, respectively, 
with d a variational parameter. For the 
repulsive potential, HD and SD, f(r) 
corresponds to the ground-state two-body 
wave function. On the country, in the case of 
the PP interaction, f(r) is chosen as the 
two-body wave function corresponding to 
the lowest excited state with positive energy. 
The pseudopotential V pp imposes the 
following boundary condition on the pah- 
wave function for zero inter particle 
distance: 



[f(r)-In(qr)ef'] r=o In{qa 2D ) 



(12) 



The free two-particle problem in 
2D-H 2 /mV 2 f(r) = h 2 k 2 /mf(r), 
holding r ^ 0 gives the general solution 



(10) F(r)=AJ„(kr)+BY () (kr) 



(13) 



in terms of the J 0 (x) and Y 0 (x) Bessel 
functions. The boundary condition (12) is 
imposed by the relation 



a 2D 



2e r 

exp 



2B y 



(14) 



In addition we require that f(r) possess only 
one node for distance 0<r<d. The boundary 
conditions at r = d together with Eq. (14) fix 
the values of the constant A and B and of the 
wave vector k. if d»a 2D , the eigen energy 

h 2 k 2 1 m is small and the node of f(r) is 
located at r =a 2D . 
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3. DISCUSSION OF RESULTS 

In this paper, we have presented 
method of evaluation of energy per particle 

(E/N) as a function of gas parameter Na\ D 

The calculations have been 
performed for three model potential V HP ,V SD 
and V pp . The resulats and shown in table Ti. 
The calculations are based on the theoretical 



formalism of Mazzanti et.al. using the 
hyperfine head chain approach and VMC at 
high densities by Xing 17 using theVMC 
method. The DMC results for the HD system 
are compared with the mean-field results eq. 

(1). From the evaluated results we observe 
that the agreement upto large value of Na 2D 
is quite good. 



Table Tl 

Energy per particle E/N corresponding to the model potential V HD (DMC results),V SD (DMC 

results) and V pp (VMC results) 

[energies are in units of h 1 1 2Ma 2D ) 



Na 2 2D 


HD 


SD 


PP 


lxlO" 7 


7.738(3) x 1CT 8 


7.732(6) x 10" s 




3xl0 -7 


2.500(1) x 10" 7 


2.498(1) x 10" 7 




lxlO" 6 


9.103(5) x 10" 7 


0.909(1) x 10 -6 


9.162(7) x 10" 7 


3xl0 -6 


2.982(2) x 10" 6 


2.997(4) x 10" 6 




lxlO" 5 


1.108(2) x 10" 5 


1.106(1) x 10" 5 


1.115(1) x 10" 5 


3xl0 -5 


3.710(4) x 10" 7 


3.702(2) x 10" 5 




lxlO" 4 


1.417(1) x 10" 4 


1.415(1) x 10" 4 


1.428(1) x 10" 4 


3xl0 -4 


4.922(5) x 10" 4 


4.910(2) x 10" 4 




lxlO" 3 


1.981(4) x 10" 3 


1.972(1) x 10" 3 


1.994(1) x 10" 3 


3xl0 -3 


7.340(2) x 10" 3 


7.724(3) x 10" 3 




lxlO" 2 


3.296(5) x 10" 2 


3.057(2) x 10 -2 


3.304(4) x 10" 2 


3xl0 -2 


1.436(2) x 10" 1 


1.105(1) x 10" 1 


1.369(6) x 10" 1 


lxlO" 1 


8.940(2) x 10" 1 


4.161(4) x 10" 1 





CONCLUSION 

We have carried out a study of the 
ground-state properties of a 2D 
homogeneous Bose gas with quantum Monte 
Carlo techniques. The universal behavior of 
the equation of state for small values of the 



gas parameter has been investigated by 
analyzing both mean-field and beyond 
mean-field contributions. By using a zero- 
range pseudopotential to describe 
interatomic interactions we have estimated 

na\ D = 0.04 is the critical density at which 
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the system becomes unstable against cluster 
formation. These results are relevant for 
present and future experiments with 
ultracold Bose gases is highly anistropic 
harmonic traps. 

REFERENCES 

1 . Proceedings of Euroschool on Quantum 
Gases in Lower Dimensions, Les 
Houches, edited by L. Pricoupenko, H. 
Perrin and M. Olshanii, J. Phys. IV. 
(2003). 

2. P.C. Hohenberg, Phys. Rev. 158, 383 
(1967). 

3. V.L. Berezinskii, Zh. Eksp. Teor, Fiz. 
59, 907 (1970) [Sov. Phys. JETP 32, 493 
(1971)]; J.M. Kosterlitz and D.J. 
Thousless, J. Phys. C 6, 1181 (1973). 

4. D.J. Bishop and J.D. Reppy, Phys. Rev. 
Lett. (PRL) 40, 1727 (1978). 

5. M. Schick, Phys. Rev. A 3, 1067 (1971). 

6. A. I. Safonov, D .S. Petrov, B. P. Vanzyl 
and N. L. Smith., Phys. Rev. Lett. (PRL) 
81,4545 (1998). 

7. A. Gorlitz, S. Berger, C. H. Schuneck, 
Y. Shin and W. Ketterle., Phys. Rev. 
Lett. (PRL) 87, 130402 (2001) 

8. D. Rychtarik, B. Emgeser, H.C. Nagerl 
and R. Grimn, Phys. Rev. Lett. (PRL) 92, 
173003(2004). 

9. N.L. Smith, J. D. Reddy, M. E. Gehm, J. 
E. Thomas and D. S. Jin, e-print cond- 
mat/04 10101 (2005). 

10. V.N. Popov in Functional Integrals in 
Quantum Field Theory and Statistical 



Physics, Reidel, Dordreht, (1983). 

11. D.S. Petrov, M. Holzmann and G.V. 
Shlyapnikov, Phys. Rev. Lett. (PRL) 84, 
2551 (2000). 

12. C. Gies, B.P. van Zyl, S.A. Morgan and 

D. A.W. Hutchinson, Phys. Rev. A 69, 
023616 (2004); L. Pricoupenko, Phys. 
Rev. A 70, 013601 (2004). 

13. D.S. Petrov and G.V. Shlyapnikov, 
Phys. Rev. A 64, 012706 (2001). 

14. T. Lofros, M.Veillete, E. G. Moon, J. 
Carbon and S. Reddy., Phys. Rev. Lett. 
(PRL) 88, 173201 (2002). 

15. M. Olshanii and L. Pricoupenko, Phys. 
Rev. Lett. (PRL) 88, 010402 (2002). 

16. F. Mazzanti, A. Polls and A. Fabrocini, 
e-print cond-mat/04 10690. (2002) 

17. Lei Xing, Phys. Rev. B. 42, 8426 (1990). 

18. S. Giorgini, J. Boronat and J. Casulleras, 
Phys. Rev. A 60, 5129 (1999). 

19. D.F. Hines, N.E. Frankel and D.J. 
Mitchell, Phys. Rev. A 68, 12 (1978); 

E. B. Kolomeisky and J.P. Straley, Phys. 
Rev. B 46, 11749 (1992); A. A. 
Ovchinnikov, /. Phys.: Condens. Matter 
5, 8665 (1993); J.O. Anderson, Eur. 
Phys. J. B. 28, 389 (2002). 

20. J. Boronat and J. Casulleras, Phys. Rev. 
B 49, 8920 (1994). 

21. R. Combsescot, A. Recati, G. Lobo, 

F. Chevy and M. M. Parish., Phys. Rev. 
Lett. (PRL) 100,180402 (2007). 

22. P. Massighan, G. M. Run, H.T. F. Stoof, 
W. Ketterle and R. G. Hublet, Phys. Rev. 
A 78,031602 (2008). 



Journal of Pure Applied and Industrial Physics Vol.2, Issue 4, 1 October, 2012, Pages (403-484) 



